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SUMMARY 

It has been shown that disturbances in external flow can significantly affect, by as 
much as an order of magnitude, the turbulent mixing rate in free shear layers. A partic- 
ularly important finding is the fact that the length scale of the external flow disturbances 
is as important as the amplitude. Also a single parameter correlates the change in 
entrainment rate remarkably well. The difference between the effect of wide -band and 
narrow-band disturbances is stressed. The inclusion of the model for pressure fluctua- 
tion term in the kinetic energy equation in a two-equation model predicts the reduced 
spreading rate in high Mach number, high Reynolds number, adiabatic, free turbulent 
shear layers. 


INTRODUCTION 

The free turbulent shear layer is a relatively simple flow which is useful for veri- 
fication of turbulence modeling and has many practical applications such as design of jet 
engine combustors, slot injection, and gas -dynamic lasers. The requirements of these 
devices generally vary. The combustor and gas laser require fast mixing and slot 
flows need slow mixing to be effective as thermal shields or to provide drag reduction. 

Of importance in engineering design of such devices is the determination of the 
relative importance of various parameters upon mixing rate, and whether one can "con- 
trol" the entrainment. The rate of free turbulent mixing primarily depends on the level 
and scale of turbulence in the shear layer. There are many things that can affect the 
structure of free turbulence; for example, large rates of strain and impinging shocks, etc. 
The purpose of the present investigation is to determine the effect of free -stream distur- 
bances when both free -stream turbulence intensity and scale are varied, and also to 
determine the influence of Mach number on the entrainment rate of free shear layers. 

The influence of free -stream disturbances may be particularly important in groimd sim- 
ulation of combustor flows, especially scramjets, where the free -stream flow is heated 
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by a variety of devices such as arc or vitiated burners which probably produce fairly ’ 
large disturbance 'revels and scales. (See fig. 1.) > . . • . . 

Ebq)erimentalists often fail to measure the intensity of such free-stream' distur- 
bances, and very rarely measure the dominant or characteristic scale. High Mach number 
is also found to be a very important factor that affects turbulent free mixing structure. 

The 1972 Langley Conference on Free Turbulent Shear Flows (ref. 1) indicated that the 
data on spreading rate for turbulent free shear layers contain disagreement between data 
and "theory" of up to a factor of 3 or more for supersonic Mach numbers, the larger dis- 
agreement occurring at the higher Mach numbers (M = 3 to 5). 

External or free-stream disturbances can be roughly divided into two important 
classes. The first class is generally referred to as acoustic disturbances. These are 
isentropic pressure waves having the following characteristics; The magnitudes of pres- 
sure, density, and velocity fluctuations are related to each other isentropically (up to mod- 
erate amplitude). These disturbances are highly directional, propagate with the speed of 
soimd, and can be transmitted and reflected by the shear layer. This type of disturbance 
is known to affect laminar shear layer and combustion stability. 

Extensive reviews of existing literature on acoustic disturbances can be found in 
Rockwell (ref. 2) and Borisov and Rozenfel'd (ref. 3). Important experimental data 
(ref. 4) showing the effect of narrow -band acoustic mode disturbances on fully turbulent 
shear layers are shown in figure 2 and are discussed in detail in the main body of the 
paper, particularly in reference to turbulence control. 

The second important class of external disturbances is usually referred to as 
free-stream vorticity or "turbulence." Both of these classes of disturbances can be 
either wide band or narrow band. Wide -band disturbances generally have typical turbu- 
lence spectra. Vinogradov et al. (ref. 5) found the core length of two-dimensional co- 
flowing jets to be strongly affected (factor of 4 difference in spreading rate) by the type of 
screen applied upstream of the external flow. Rodi (ref. 6) also reports similar findings 
(Vagt (1970) and Patel (1970)). The present paper describes a study of the effect of wide- 
band disturbances, employing an equation for the scale of turbulence since disturbances in 
the "free-stream" (on either or both sides of the mixing layer) would generally have disr ^ 
tinct scales which are not related to the local scale of the shear layer; A parametric 
study of the effect of scale and intensity of free-stream disturbances is presented so that, 
for example, one can assess the influence of tunnel disturbances on free -flight simulation 
of scramjet engines. (See fig. 1.) 

The present numerical solution of the free turbulent shear layer employs a "two- 
equation" turbulence model, a turbulent kinetic energy equation for the intensity, and a 
dissipation equation for the length scale of turbulence. A model for the pressure -velocity' 
correlation term, representing a "compressibility" effect, which was developed by Oh 
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(ref. 7), based on an eddy -shock -wave concept is included in the turbulent kinetic energy : 
model equation for predictions of the influence of high Mach number upon entrainment in 
free shear layers. 

SYMBOLS 


a 

Cp 


mass fraction of species 
specific heat 


e , turbulent kinetic energy, -i(u"2 + v"^ + w"2) 

£t 

f frequency 


H 

L 

I 

^in 

M 

n 

^Pr,e 


length scale 
integral length scale 

Mach number; also grid size 

total number of species; also exponent of isotropic turbulence energy initial 
decay law 

ratio of turbulent diffusivity in turbulent kinematic energy equation to that of 
mean momentum 


n 

total enthalpy, ^ 


i=l 


reference length 


Npr,t 

Npr,e 

P 

R 


turbulent Prandtl number 

ratio of turbulent diffusivity in dissipation equation^to that of mean momentum 
static pressure 

gas constant; also Reynolds number 
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Nsc,t 

T 

u 


Us 


V 


w 



X 


Xc 


Xo 


y 

y 5 

Y 

6 

6o 


e 

A 


I 

i 

f 

I 

i 

tiu'biilent Schmidt number 

I 

j 

i 

static temperature 

velocity component in x -direction 

I 

1 

local sonic speed 
velocity component in y-direction 
velocity component normal to u and v 
virtual origin of velocity similarity profiles 

i, 

coordinate parallel to outside flows (fig. 1) 

I 

core length of jet 

i 

1 

initial x where initial conditions are specified 
coordinate normal to x (fig. 3) 

r 

y where (u - U 2 )/(ui - ^ 2 ) = 0.5 

= (m 2 - l)^^yM 

ratio of specific heat 

shear layer width parameter, ly^-o 1 - Yu=o 9 ]/^ 


= 6 at X = Xr 


dissipation rate of turbulent kinetic energy 


computed shear layer width referred to that implied by 
with computed a and Xq, Yu=/ 0 ' yS=\/6r9|/[^- 


/ 


linear spreading law 
32(x - Xo)/aj 
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\ = 0.8756o 

fie eddy viscosity 

p density 

(T similarity parameter used as a measure of spreading rate 

I 

I 

oq low-speed constant -density flows 

CT* = lim a 

/e^oo - 0 

<p = y . 

Ui - U2 Uj - U2 jp.QL^Jx=0 

Subscripts and superscripts: 

1 conditions on high-velocity side external flow 

2 conditions on low-velocity side external flow 

i species i 

i,j tensor indices 

max maximum value 

°° ' conditions outside of shear layer 

(“) conventional temporal mean 

( )’ ' (instantaneous) - C”) 

(~) . mass -averaged temporal mean 

’ 1 ■ 

( )" (instantaneous) - ('“) 
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GOVERNING EQUATIONS AND ASSUMPTIONS 


The physical problem is depicted schematically in figure 3. The shear layer is 
formed between two parallel uniform flows of different velocities. The higher velocity 
is called u^ and the slower stream velocity is denoted by U 2 . The X-axis is alined 
with the external imiform flows and the Y-axis is the direction normal to the uniform 
flow. (The flow properties vary most rapidly in the y-direction,) In the subsequent 
^analysis, mass-averaged dependent variables, proposed by Fayre (ref. 8), are used for 
most flow properties. Mass -averaged values are represented by a superscript tilde ' 
("'). The following relationships hold: . ■ 

~ _ p'u’ 

u = u + _ ■ 

p . 



so that ^ Both independent and dependent variables are used in their primitive 

'form. The governing equations employed cover both multispecies and high-speed 
mixing. 


Simplifyir^ Assumptions 
The following assumptions are made: 

(1) The flow is steady on the average 

(2) The Reynolds number is sufficiently high so that the mixing layer is fully tur- 

bulent and molecular diffusion can be neglected compared with turbulent 
diffusion 

(3) The axial pressure gradient is zero 

(4) No chemical reactions occur 


Basic Equations 

The typical free shear layer with no extraneous strain has an, order of magnitude 
difference in the dierivatives of the' flow properties in the x- and y-direction and allows the 
usual boundary-layer -type (quasi -parallel) approx imationsi The conservation equations 
describing the problem for the mean flow quantities may be written in the following form 
(ref, 8): ; ' ^ 

Continuity: 

^(y^) - ff ' (j) 

dX^.dy 
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”x" mean momentum: 


dK dy ’ by 


•y" mean momentum: 


8p 9(pv”2) _ 

by by 


Mean energy: 


pu -^ + pv - y~J 

9x ®y ^ 3y 


3(yJpBTv^) _ Q 


Species: 


9as ^ 3a j T 3(yJpaiV’) 
pu — i + pv — i + y-*' ^ — i 

3x 3y 3y 


Equation of state: 


P = 




( 2 ) 


( 3 ) 


( 4 ) 


(i = l, 2, . . n) 


(5) 


( 6 ) 


(yi-l)f^^^^Riaij 

In these equations all the variables are rendered dimensionless by referring all lengths 
to L, velocities to uj, pressure to PiUj^, density to pj, temperature to ui^/Cp j, 

total enthalpy to uj^, specific heat to Cp^i jcp^l = ^ Cp,l,iaijij, and gas constant to 
1^1 ^1 ~ ^ ^l,i^l,ij* The exponent J»0 for two-dimensional and J » 1 for 

\ i=l ’ 7 

axisymmetric shear layers. 


Closure Assumptions 

The preceding set of equations are not closed because the turbulent correlations 
are unknown. Boussinesq's eddy viscosity model is assumed to describe the turbulent 
shear stress; that is, pu"v" = -\i^^jby along with Prandtl's energy model for the eddy 
viscosity, pg = Culpe, where Cp is a modeling constant and e is the turbulence 
kinetic energy e = p(u"2 + v"2 + w"2)y^. Also pH”v"/(bS/by) and p^(v" J {b^i/by) 
are assumed to be proportional to the eddy viscosity (by constant factors); that is. 
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pa’V' - 


where 


_ f^e 3a _ /^e „ af , *^e Y , 7 ^,, 9u 

^Sc,t Npr,t turbulent Schmidt and Prandtl number, respectively. 

Turbulence Model 


Prandtl' s energy model for the eddy viscosity requires the turbulence kinetic 
energy and length scale. Often the length scale is algebraically related to the width of 
shear layer itself and yields a reasonable result for the free shear layer. But, in order 
to study the effect of external disturbances which have their own scale (which is not 
related to the shear layer), use of an equation governing evolution of the scale of turbu- 
lence is essential. Thus a "two -equation" model approach is the minimum necessary to 
describe the turbulence field in the present problem. The equation for the turbulence 
dissipation rate is used to compute the turbulence scale. There are other equations 
which would serve the same purpose, such as a vorti city equation (ref. 9) or an equation 
for the quantity, energy -times-length scale (ref. 10). Actually, these equations are all 
very similar (ref. 11 ). 


Equations derived from the Navier-Stokes equations to describe turbulence energy 
and dissipation include a number of unknown correlations. Modeling of these unknowns 
requires a delicate balance between mathematical rigor and physical intuition since elab- 
orate mathematical manipulation can be meaningless if not supported by experimental 
measurement. In general, the higher order correlations in the turbulent kinetic equation 
are better known than those in the equations for the length scale (ref. 12). Most notori- 
ously difficult to measure, and hence least known of the turbulent fluctuating properties, 
are the pressure fluctuation terms. A model for the pressure -velocity correlation 
hypothesized by Oh (ref. 7) is included in the kinetic energy equation. No attempt was 
made to model pressure fluctuation terms in the dissipation equation. The model used 
by Spalding's group with their constants (ref. 1) was used herein for the dissipation equa- 
tion. The final model equations are as follows: 

Turbulent kinetic energy: 


pu — + pv — - 
ax ay. 



^Pr,e ^57 



+ Ce^ ipOe + p€ 


- Ce,2H(M - me 


^tslgn®U8^|^-0uil||| 


(7) 


Turbulent dissipation equation: 


^ ^ 91 - il\ _ C. 1 1 = 0 

ax ^ ay y ay^y^ ay) "'e.l g '^®l9y; * q 


(8) 
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where J = 0 for two-dimensional and J = 1 for axisymmetric shear layers, H is the 
Heaviside unit fimction H(0 = 0 when ^ < 0 and H(?) = 1 when ? ^ 0. 'Also, the 
following definitions are used: D = aii/0x + 3v/9y; Us = Local sonic speed; 

/3 s (m^ - l) so that Pe = The constants used are 

Cfi - 0.09, a 0.7, a 0.9, - 1-0, Hpj,^g = 1.3, Cg^l * 0, 

^e,2 “ 0*14, Cf 1 a i.43j and Cg 2 = 1-92. 

Initial and Boundary Conditions 

Initial conditions , r Initial profiles of all dependable variables are required to start 
the solution procedure. All the results reported in this paper are computed with a set of 
initial profiles which were "guessed at" based iq)on the physics of the flow. The scarcity 
of experimental data, especially measured profiles for the turbulence energy and length 
scale, is the reason for using "arbitrary" (but physically reasonable) Initial profiles. 

The quasi -similarity (all the computed results show linear shear layer growth, see fig. 5) 
of the results is a posteriori justification, that is, the calculations proceeded far enough 
downstream for the influence of the starting profiles to be "washed out." The initial pro- 
files are generated as follows: First axial velocity u profiles are generated for a 

given width parameter 6 q = ■■ ^ — fl-0*9|x=0 

L 

u(0,y) = ij(ui + U2) + (ui - U2)erf(?)J (9) 

where ? = 2 x 0.90621514y/6o. Then 



®fr ®sh 


where the first term e^j. is the contribution of free -stream disturbances and the sec- 1 
ond term is chosen as the contribution of mean shear; ^rom the assumed relation 

e « when gj » §2 * The term ejnax,sh another input which is chosen as 

0.12 for most incompressible cases and as 0.08 for s^ersonic cases. The term ig}^ 
was neglected when e^^ -f 62 > 2ejQ3x,sh* these equations are in nondimens ional 

form. In order to estimate dissipation, the distrilxition of length scale is computed as 
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( 11 ) 


2(0 


^ ^[(^1 ^2) + 1^1 - i2|®^^(?)j^r >.emax,sh - ^(Zi + Z2)(®1 + ®2) 6gjj 

® ®max,sh - 4(®i + ®2) ^ 




2 fr 


^sh 


where A. s 0.8756o and 2gh = 0 when emax,sh ^(©1 + e^. Then e is computed as 
6 = 53 / 2 / 2 . Total enthalpy is assumed to be initially uniform 

H(0,y) = l 

The static temperature is obtained from the expression 


f (0,y) = 




m.2ih _ y -.,1 M^2{i2 


Cp(y - l)R 


(12) 


For static pressure, normal velocity and density, the initial assumptions include: 
P(O.y) = Poo " 
v(0,y) = 0 i 

p(0,y) = Mj2 ^ 

TR J 


(13) 


Boundary conditions . - The external free stream usually contains some type of dis- 
turbance. The major disturbances can be divided into two broad groiq>s, acoustic and 
vortical. The effect of narrow-band acoustic disturbances is fundamentally different than 
that of free -stream vorticity (turbulence). The present research effort is concerned with 
the effect of wide -band -type free -stream turbulence only. 

The free -stream turbulence could be the result of upstream agitation, diffusion 
from an acyoining shear layer (wall boundary layer in enclosed flow or free shear layer 
in case of coaxial jet, etc.), or distributed sources such as chemical reaction or distrib- 
uted fine obstacles, etc. (that is, dust flows). The most well dociimented data are for 
decaying Isotropic turbulence. (A good example is grid -generated turbulence.) There- 
fore, the boundary conditions for the kinetic energy and dissipation equation are limited 
to isotrqpic decaying turbulence in the present study. 

Batchelor and Townsend (ref. 13) found grid -generated isotropic turbulence decays 
in the early stage as 
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^in _ c(^ 

M -- \M ■ M/ 


1/2 


where A, B, and C are constants, and Ijn, Xg> and M are the integral length scale, 
grid location, and grid size, respectively. ' 

In order to study the variation of boundary conditions for the turbulence model 
equations with x, the turbulence energy is-assiuped to vary as ' 


e«, = UoO^Al:^ - ^ 

M M 


n 


In the absence of mean shear, the model equation for the turbulent kinetic energy equa- 
tion becomes 


dx 


Substituting the expression for Coo results in 

A 

t = 5 3 JL A iL 


From the relation I = 


_eV2 




/_ MAV^ fx ' ^gV 

n \M ■ M/ 

Initial boundary conditions are therefore 




n 




n \ = MA V^/ 'xq fg' 

«>/x=0 ■ n Im ■ M 


From equations (14), (15), (17), and (18), 
®oo = (Ojj=0*^ ■ 


1 

■ X - Xq ■ 

n 

^®~^x=0 




u 2 

J 


(14) 


(15) 


(16) 


(17) 


(18) 


(19) 
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( 20 ) 


^OO “ 


fen.. 


f- 

X - Xq 

n-1 

(®«’^x=0 

n-|il 

Wx=.o" 


1 

r 

to 

i . 

1 ^ 

/ 


The recent measurements made by Comte-Bellot and Corrisin (ref. 14) show that 
n « -1.25. 


Symmetric boundary conditions are imposed on the X-axis in the axisymmetric and 
two-dimensional jet cases. (Center velocity Is allowed to decrease.) 

Nunierical solution procedure . - The numerical method used tp solve the governing 
equation is the implicit finite -difference method of Crank -Nicholson with the dependent 
variables in conservation form. Linearization is accomplished by initially lagging the 
nonlinear coefficient a step and then integrating until convergence. 

• Variable grid spacing is used in the y-direction. The spacing was increased by a 
geometrical progression on each side of the shear layer center line. In order to prevent 
the shear region from outgrowing the computing net, the computing net was continuously 
expanded in the y-directiqn by doubling the grid spacing whenever the converged solution 
for u changes by a predetermined margin at both edges for the pure shear layer case 
and at the outer edge for the jet case. A t37pical rim, 250 cross-node points and 
170 marching steps in x («100 times the initial thickness distance), for the incom- 
pressible binary mixing case takes approximately 200 seconds of CPU, time on the 
CDC 6600 computer system. Sixty percent of this time is actually spent in solving the 
finite -difference equations and the rest of the machine time is spent in data management 
and plotting routines. A solution for the supersonic shear layer takes approximately 
three times longer than for an incompressible case. Convergence was tested on u, v, 
and e with a convergence criterion 0.1 percent relative (or combined with an equivalent 
absolute criterion for v). 

RESULTS AND DISCUSSION 
Comparison With Data 

Experimental data which can be used tp check the validity of the present prediction 
are extremely scarce. Most often, data are incomplete, that is, the intensity of free- 
stream disturbances is measured but not the characteristic scale. Also, the core length 
of the jet may be measured but not the initial boundary flayer thickness (ref. 5), etc. 

Thus the comparison with incomplete data, becomes no more than finding the missing 
data which would 3 ne Id agreement, in the results. Nevertheless Rodi (ref., 6) demonstrated, 
by using a prediction method similar to the present one, that the measured effect of free- 
stream disturbances (Patel (1970) and Vagt (1970)) on free turbulent mixing could be qual- 
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itatively simulated. In this vein, a comparison will be made with the Vinogradov et al. 
(ref. 5) ejq)eriment. The experimental' setup (see fig. 4) consisted of an enclosed two- 
dimensional channel with two separating plates which initially divided the channel into 
three equal size ducts. The ratio of center duct velocity to the velocity of two outer 
ducts was varied along with the "agitating grids" in the outlet sections of the ducts. The ' 
experiment determined the variation of core length of the center jet with velocity ratio 
,^d initial. intensity of turbulence for incompressible isothermal flow. Frorn the present 
authors' viewpoint, the data have the following defects: (1) The wall shear layers on eithei 
side of the splitter plate were not measured. (The core length depends not only on the 
spreading rate but also on the virtual origin, that is, initial shear layer thickness.) 

(2) The core len^h was measured from schlieren photographs, with smoke in the center 
jet. . However, the end of the core was not defined clearly, that is, whether the concen- 
tration was ax = 0.99 or 0.9, etc. The analx^ical determination of Xc depends greatly 
on such a criterion. (3) The location where the grids are placed is not certain. If the 
grids were placed at the exit of the jet as the article (ref. 5) implies, then the results 
would have been considerably affected by the immediate wake of the grid. (4) The outer 
uniform flow carries not only decaying grid -generated turbulence but also turbulence - 
diffused from the outer wall boundary layer. (The nature of such turbulence is much 
less well known than grid turbulence.) (5) The initial free -stream turbulence energy 
was measured but the scales were not. However, because of the lack of more complete 
data, comparisons will be shown for this case by using plausible assumptions. 

A number of predictions were made with arbitrary but reasonable initial shear 
layer widths for the no grid case. (See circle S 3 onbol in fig. 4.) The prediction and 
experimental data are compared in figure 4. The flow quantities assumed in predictions 
of all cases are ui = 30 m/sec, 6q = 0.5875 x 2 cm, definition for the core length 
ai = 0.98, length scale for the no grid case fooA = is assumed to be pro- 

portional- to the grid spacing (or rod diameter). The implication of this comparison is 
the slc^e of the no grid case data with velocity ratio is used to define the end of core as 
aj = 0:98, the Xg at ^2/^1 = 0.3 for the no grid case is used to calibrate the initial 
shear layer width 6 q = 0.5875 x 2, and the grids with 2-mm rod and 4 -mm rod are used 
to fix A(iooA)/(2-mm rod diam) = 0.15. Therefore, the results for the 8-mm grid can be 
looked upon as the real comparison, which is certainly very reasonable and tends to vin- 
dicate the present approach. ' ' 

Parametric Study of Effect of Free-Stream Disturbances 

Because of the lack of complete data, the comparison of the present results with 
experiment is quantitatively inconclusive (but qualitatively satisfactory) as the attempt 
in the previous section demonstrates. ' 
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In this section a parametric study is made of the effect of the intensity and length 
scale of wide -band -type free -stream disturbances on 2-D incompressible, free, turbulent 
air -air mixing, A number of computer rims were made with identical flow conditions 
except for variations in velocity ratio ^ 2 /^h turbulent intensity \/e^/uoo, and length 
scale "Zoo/X.. All computed cases are summarized in table I with the results a and Xq. 

To show the effects of Coo and on profiles of shear layer properties, 
samples of the similarity profiles u, e, and l/\ as functions of o’y/(x - Xq) and 
^max> ®max> ^max> and A as a function of x (which represent large disturb- 
ances, very small disturbances, and effects of disturbance on both sides of the shear 
layer and on just one side alone) are shown in figure 5. The spreading parameter o 
which is used as the measure for the spreading (see fig. 3) is defined herein by the rela- 
tion (ref. 15) a « 1.32/Atj where Atj is the angular distance between two rays when 
(u - U2)/(Si - U 2 ) =* (0.1)^/^ and (u - U 2 )^(ui - U 2 ) = (0.9)^/^. The accompanying virtual 
origin is Xq. The quantities a and Xq are evaluated from the last u profile com- 
puted and a profile at about X » 0,8 of the total distance computed. A single parameter 
which can be used to compare the rate of spreading of all properties in all flow conditions 
is not yet known. The par;ameter a and other parameters (ref. 16) have their lim- 
itations. The quantity a ccan be used to compare the growth rate of velocity width (also 
momentum thickness and entrainment rate for constant -density shear layers) only when the 
velocity profiles are similar. But the presence of a high level of turbulence in the external 
flow distorts the profiles; however, they are "self -similar.” The self-similar velocity pro- 
„ file for the quiet boundary condition case (fig. 5(a)) is similar to that of Liepmann and 
Laufer (ref. 17). (A comparison can be found in ref. 7.) The presence of high ej makes 
the self -similar velocity profile smoother near the high-velocity side because of the 
increased shear stress. (See fig. 5(e).) The presence of high free-stream disturbances on 
both sides stretches the profile at both ends so that the self-similar profile becomes nearly 
similar again. (Compare figs. 5(a) and 5(i).) The a values used in table I, figure 6, 
and figure 7 can be regarded as comparing approximately 80 percent of the momentum 
: thickness because of the way it is defined. 

In the course of this study, the following qualitative observations were made: ' 

(1) The effects on turbulence intensity is approximately proportional to \/e^y/(ui - S 2 ); 
that is, for a given \/e^Uoo, the effects of disturbances in the high-velocity side became 
more prominent and the effects of external turbulence on the low -velocity side diminishes 
as U 2 /uj decreases. The higher U 2 /ui, the more sensitive the flpw is to external dis- 
turbances, (2) As Zoo -- 0 \f^ — 0, ^ and Xq approach an asymptotic value. 

(3) The overall profiles are quasi-similar (i.e., though turbulence properties are still 
changing the mean velocity profiles are similar) in most of the cases computed even 
though the boundary conditions imposed are strongly nonsimilar. Note the constancy of 
^ “ |yu=\/0 “ yfl=\/079|/[^-32(* " ^o)/^ ^ function of X (linear spreading) in fig- 
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ures 5(d), 5(h), 5(1), and 5(p). (4) Very small causes numerical instability near the 

‘ ■ s3/2 ■ 

edge due to the definition I = (5) The ratios of diffusivities in e and e equations 

have to have a certain value, that is, 1.3 (as Spalding’s groxq) used). Otherwise the lei^h 
scale profile would exhibit an anomalous dip or peak near the edge of the shear layer. 

(6) Disturbances on the low -velocity side of the external flow can further increase mixing 
rate up to 10 percent over the correlated results shown in figure 6. 

Surprisingly, a single parameter 

^ _ ui + U2/yg^ j /ZA 

ui - a2\ui - U2/x=o ^^x=0. • 

correlates the variation of o/o^ with "one side" external disturbances as shown in 
figure 6. . ' . 


Here a* is the limiting value of a for small (p: (Actual a* values used are 
marked in the table.) It should be noted that the present calculations, as correlated in 
figure 6, indicate potentially large effects of free-stream disturbance on spreading rate 
for simple shear layers (factor of up to 10 change in spreading rate possible). 






/ i/ei ^1) 

In figure 7, predictions made with — :: — -I = 1.7 and 0 are plotted with 

W-"2Vx=0 

experimental data collected by Birch and Eggers (ref. 1) in the format suggested by Kline 
for <T as a function of U2y/ui. This comparison suggests that the disturbances in the 
facility free stream were responsible for at least some of the data scatter which is more 
than 100 percent at the high velocity ratios ^U2/ui 


The effect of narrow -band acoustic disturbances could be markedly different from 
the effects of the wide -band vorticity type so far considered. The difference is very well 
demonstrated in the experiment conducted by Vlasov and Ginevsky (ref. 4). The experi- 
mental setup and the results were shown in figure 2. Narrow -band acoustic disturbances 
generated by a loud speaker were injected into a fully turbulent free subsonic jet. The 
center-line velocity decay with x is plotted for the no-sound case and two different dis- 
turbance frequencies. The data show that the external disturbances not only accelerate, 
mixing (low-frequency disturbance), but could also attenuate mixing (high-frequency input). 
This is a very significant result and indicates that the mixing rate can be controlled by 
artificiaUy inputting narrow -band disturbances. It is currently planned to study this 
phenomenon with a spectral plane analysis. 


■*^Note that the slow rate of development to a "pure shear layer" is also conjectured 
(ref. 16) for the large scatter of data as U2/^l 
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Influence of Mach Number 

Oh (ref. 7) modeled the pressure -velocity correlation i term p' in the kinetic 

energy equation based on an eddy shock -wave concept. The prediction (ref. 7) made by a 
one -equation model was reasonably good. The same cases were computed again with the 
dissipation equation included. (Note that the p* terms are not yet included in the dissi- 
pation equation.) 

A prediction with Cg 2 = 0-34 (equivalent to C 5 C 0 = 0.07 of ref. 7) are shown in 
figure 8 .' The plots of a against Mi are qualitatively similar to the results of the 
one-equation method (ref. 7). Detailed profile comparisons indicate the need for slightly 
higher diffusivity in the dissipation equation which is a possible indication of the need for 
a model of the p' term in that equation. 

To be noted is that all the data shown in figure 9 are from relatively "clean" con- 
figurations. The data that have some sort of wall proximity, which may allow reflec- 
tion of large noise levels back into the shear layer or allow self generation of low- 
frequency narrow -band noise such as cavity flows, are excluded. These "noisy flows" 
generally spread faster. (See fig. 9.) 

CONCLUDING REMARKS 

It has been shown that disturbances in external flow can significantly affect, by as 
much as an order of magnitude, the free turbulent mixing rate in shear layers. A parti- 
cularly important finding is the fact that the length scale of the external flow disturbances 
is as important as the amplitude. 

Also, a single parameter correlates the change in entrainment rate remarkably well. 
The difference between the effect of wide -band and narrow -band disturbances is stressed. 
The inclusion of the model for the velocity-pressure correlation term in the kinetic energy 
equation in a two-equation model predicts the reduced spreading rate in high Mach num- 
ber, high Reynolds number, adiabatic, free turbulent shear layers. 


y 
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TABLE I. - COMPUTED RESULTS OF SPREADING RATE AS A FUNCTION 
OF FREE-STREAM DISTURBANCES , , 


Plow conditions for data given in table, 
uj = 30 cm/sec 

Ttotal = 318 K , 

Ptotol-^™ 

Ml = 0.084 
R/cm = 17296.5802 

6o = 0.3 cm = - y -,0 ’ 

■ ; ! 

All cases computed from Xq = 0 to x = 30 cm (1006 q) 
Varied conditions, at x = Xq ’ 

0.5, 0.3, 0.6, 0.9 

“1 


■J 


= 0.005 to i/0T03 


= 0.1 to 10 

A 


. (\/.= 0.8756o) 


uo \/e<» 2 ‘«> 2 / 

(a) ^ = 0.05; V, ■ >:: = 0.02; ^=0.1 


^1 


«2 


foo^l 

\ 

Value of a and for ’j- of - 

^ «l/ 

0.06 

0.086 

0.1 / 

yo.o3 

■ 1,,; 

4 

6 

10 

13.83 

-12.99 

13.69 

-12.84 

13.54 

-12.69 

ai4.14... 
-13.41 , 

13.72 
-1 5.^01 

/ 

13.10 

-l4.00 

It 

12.73 
-/l5.73 

i 

12.83 

-23.47 

12.24 

-24.23 


^Used as o* in figure 6. 




TABLE I. - COMPUTED RESULTS OF SPREADING RATE AS A FUNCTION 
OF FREE -STREAM DISTURBANCES - Continued 


Oo 2 2 

(b) ^ = 0.3; = 0.005; ^ = O.C 


Value of a and — for 

\/®oo 1 

V 'of - 


"1 


0.0377 


23.13 

-20.76 


21.92 
-21.90 

21.92 

-21.37 

21.07 

-21.03 


20.86 

-27.50 

20.05 

-27.31 

19.01 

-26.10 


19.39 

-38.55 

18.06 

-39.73 

16.26 

-38.36 


(c) ^ = 0.3 
«1 


^*”, 2 . ^°°,1 *° o ,2 


X \/®«o 1 


• 

^ . "1 


0.06 

0.1 

VOS 

V05 




^3.36 




-9.79 . 


- - : • 


22.22 


1 • • • 


-26.78 

i 

21.89 ■ 
-26.10'" 


21,05 

-33.75 

21.86 

20.74 

19.78 

19.21 

-21.98 

-27.73 

-34.40 

-39.70 

21.44 ■ 

19.92 


17.85 

-2L52 

-27.62 


-41.03 

21.01 

18.88 

17.10 

15.99 '■ 

-21.26 

-26.63 

-34.13 

-39.62 


^Used as a* in figure 6. 


























TABLE I. - COMPUTED RESULTS OF SPREADING RATE AS A FUNCTION 
OF FREE-STREAM DISTURBANCES - Continued 


(d) 




0.005; 


too 2 


= 0.1 



^Similarity profiles are included in figure 5. 



.(e) 

=2 = 0.6; 

V^»,l _ 

6oo 2 

1 _ ^*.2 




“2 ’ ^ 

■ '■ ^ 


B 


Value of 

a and ^ for 

oo i 

of - 

>1 


H 

IQQI 

0.6 


/oloTs 

vC02 

vC03 

0.1 

BI3 

; .. 

1 ■ * 




1 

_ 

. 

. 

43.86 





• 


-54.36 



4 


39.44 

i 34.00 



28.52 



-45.60 

-54.59 



-87.43 

10 

42.20 




*22.79 

*19.71 , 


-43.92 




-63.58 

.98.48 


^Similarity profiles are included in flgpire 5. 









































TABLE I. - COMPUTED RESULTS OF SPREADING RATE AS A FUNCTION 
OF FREE -STREAM DISTURBANCES - Concluded 


(h) “2 = 0.9; = 

ui uj U2 ^ ^ 


ioo^l 

X/N v/®«> 1 

Value of a and for — *- of - 

^ Uj 

0.005 

\/ 0.03 

0.1 

209.37 



-179.20 


1 

199.09 



-168.19 


10 


23.87 



-112.38 
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SIMULATION 


Figure 1. - Schematic of Ramjet/Scramjet engine ground simulation. 
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Figure 2 .- Effect of narrow -band disturbances. 




Figure 3.- Sketch of problem. 



Figure 4.- Effect of wlde-band free -stream disturbances on turbulent free mixing 

(low speed) comparison with data. 
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(a) Variation of u/uj„ax with o(y - y 5 )/(x - Xq) for U 2 /S 1 = 0.8 with turbulence 
boundary conditions at x = 0; \/l^/Uoo = 0,005; i«o/A. = 0.1; L = ^ fio; 

Xo/L = -11.5034; 0= 49.3547. 
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(c) Variation of ai/(x -Xq) with a(y - y.5)y^(x - Xq). Xq/L = -11.5034; a - 49.3547 
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(d) Variation of Umax/3ref» emax/eref» («"v")max/^^'^")ref» ^ */^- 
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Figure 5.- Continued. 
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(l) Vairiation of fl/flniax with a(y - y 5 )/(x - Xq) for U 2 /uj = 0.6 with turbxilence - 
boundary conditions at x = 0; \/i^/Uoo = \/0.03; iooA = 10; uniform initial length 
scale profile; Xq/L = -25.8005; cr = 19.7200. 
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(j) Variation of e/emax ^(y - y. 5 )/(x - ^o)- Xq/L = -25.8005; a ■= 19,7200. 

Figure 5.- Continued. 
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(k) Variation of oZ/(x - Xq) with a(y - y.5)/(x - Xq). Xq/L = -25.8005; a = 19.7200, 
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Figure 5.- Continued. 
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Figure 5.- Concluded. 
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Figure 6.- Parametric correlation of predicted effect of wide -band free 
disturbances (only one side) on turbulent free mixing (low speed), a 
on width for 10 percent and 90 percent momentum. 



Figure 7.- Effect of wide-band free-stream disturbances. Variation of 

(ui - U2)/(ui + U2). 
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Figure 8.- Effect of Mach number on spreading rate prediction for clean flow. 
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(a) Cavity Flows (refs. 22 and 23). 



(b) Massive wall injection (ref. 24). 
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(d) Backward facing step with 
matched pressure (ref. 20). 



(e) Near field of isolated jets 
(ref. 21). 


Ml 
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M2.2 1 


(c) Triple shock intersection (ref. 5). 
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(f) Sonic slot injection into siqier- 
sonic hypersonic flow (refs. 26 
f and 27) co-axial jet (ref. 28). 

Figure 9.- Experimental configurations used for supersonic shear layer experiment. 
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